Impurity-induced polaritons in a one-dimensional chain* 
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A detailed analytical study of an impurity induced polariton band arising inside a spectral gap 
between lower and upper polariton branches is presented. Using the microcanonical method, we 
calculate the density of states and localization length of the impurity polaritons. Analytical results 
are compared with numerical simulations and excellent agreement is found. 
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I. INTRODUCTION 



An opportunity to create Anderson localized states 
of electromagnetic excitations in disordered dielectrics 
has attracted-!^ great deal of attention during the last 
two decades.Ela Though the localization of light in three- 
dimensional uniform-on-average random systems proved 
to be difficult to achieve, new opportunities opened 
up with development of the structures with periodi- 
cally modulated dielectric properties, so called photonic 
crystals.Bu It was first suggested in Ref. |^ that the pho- 
ton localization can be more easily achieved if the photon 
density of states (DOS) is depleted at certain frequency 
domains. Photonic crystals allow the creation of photon 
band-gaps, regions of frequencies where the photon DOS 
is zero. Introducing isolated defects in such structures, 
one can create local photon modes similM to. well known 
defect phonon modes in regular cazstalsErt2l or electron 
impurity states in semiconductors.tiJ A forbidden gap in 
the photonic spectrum occurs when the electromagnetic 
wave length becomes comparable to the lattice constant. 
Consequently, for microwave, infrared and visible ranges, 
the photonic crystal is a structure of macroscopic dimen- 
sions. Local states in these structures also arise due to 
defects of macroscopical dimensions. The great interest 
in photonic crystals is promoted by their highly unusual 
quantum electrodynamic properties. For instance, spon- 
taneous emission, which is completely-suppressed within 
the frequency range of the band gaps,Eican be effectively 
controlled by introducing local modes.E^lij 

The intriguing possibilities of photonic crystals initi- 
ated interest in optical effects in other types of photonic 
band-gaps. Such gaps can arise, for example, between 
different polariton branches. The fact that the frequency 
region between polariton branches can represent a stop- 
band for electromagnetic waves propagating in certain 
directions has been well known for a long time. The 
reflection coefficient in this situation can reach the mag- 



nitude of 90 percent or even greater. This effect was 
used in the 40's and 50's in order to create monochro- 
matic infrared beams (it was called the method of resid- 
ual or reststrahlen rays, and the respective spectral inter- 
val was called the reststrahlen region). In certain cases, 
however, the stop bands exist for all the directions, and 
then a genuine spectral gap arises with the same con- 
sequences for quantum electrodynamics as in the case 
of photonic crystals. This was first realized in Ref. 
where it was shown that a two-level atom can form an 
atom-field coupled states with suppressed emission simi- 
lar to the states discussed in Ref. |lj for photonic crystals. 
Periodic arrangement of such two-level systems produces 
an impurity-induced polariton band within the polariton 
gapM 

Local polariton states associated with a regular iso- 
topic defect without any intrinsic optical activity were 
introduced in Refs. |8 19. These states are coupled 
states of electromagnetic excitations with phonons (or 
excitons), with both components, including the electro- 
magnetic component, being localized in the vicinity of 
the defect. These states are, in a certain sense, analo- 
gous to local photons in photonic crystals considered in 
Refs. 13 1^, although they arise due to microscopic im- 
purities in regular crystal lattices. On the other hand, 
local polaritons can be considered as local excitations of 
a crystal coupled with the electromagnetic field. Regular 
local phonons (excitons) also interact with the electro- 
magnetic field, but this interaction results mainlyJji ab- 
sorption of light and radiative decay of the states £3 The 
local polaritons arise in the region where electromagnetic 
waves do not propagate, and there are, therefore, neither 
defect induced absorption of light nor radiative damping 
of the local states. The electromagnetic interaction leads 
in this case to new optical effects, and strongly affects 
the properties of the local states. It is well known, for 
instance, that local phonons or excitons in a 3D system 
(as well as local photons in photonic crystals) arise only 
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if the "strength" of a defect exceeds a certain threshold. 
Local polaritons in systems with an isotropic dispersion 
split off the allowed band without a threshold even in 
three dimensions £3 This effect is caused by the inter- 
action with electromagnetic field, which results in the 
Van-Hove singularity in the polariton DOS. 

One of the optical effects caused by local polaritons 
is resonant tunneling of electromagnetic wave with gap 
frequencies. This effect was first suggested in Ref. |2^, 
where the results of numerical simulation of electromag- 
netic wave propagation through a relatively short ID 
chain with a single defect were reported. It was found 
that the transmission coefficient at the frequency of the 
local polariton state increases dramatically up to the a 
value close to unity. In spite of the general understand- 
ing that local states should produce local tunneling, this 
result still seems surprising because transmission of light 
is effected by a defect with microscopic dimensions much 
smaller than the light's wavelength. In addition, the en- 
ergy of the electromagnetic component of local polaritons 
is much smaller than that of the phonon component. Tra- 
ditional wisdom based upon properties of conventional 
propagating polaritons tells that, in such a situation, 
most of the incident radiation must be simply reflected. 
The results of Ref. ^ demonstrated that this logic does 
not apply to local polaritons. The physical explanation 
of the result can be given on the basis of consideration 
of local polaritons as the result of interaction between 
electromagnetic field and local phonons. The latter have 
macroscopic dimensions comparable with light's wave- 
length, making the interaction effective. An electromag- 
netic wave is carried through a sample by phonons, which 
is tunneling resonantly due to the local state with mostly 
phonon contribution. Tunneling was confirmed in Refs. 

where we analytically calculated the transmission 
coefficient through a linear chain with a single defect. 
We found that a light impurity indeed gives rise to the 
resonance transmission through the forbidden gap of the 
polariton spectrum. The resonance occurs when the de- 
fect is placed at the center of the chain. The transmis- 
sion at the resonance becomes independent of the chain's 
length and reaches the value of unity if the resonance 
frequency coincides with the center of the polariton gap. 
The analytical calculations provided an explanation for 
the numerical results of Ref. in particular it explained 
a strongly asymmetrical frequency profile of the trans- 
mission coefficient, and the presence of a deep minimum 
following directly after the maximum. 

Analytical calculations in Ref. ^ were supplemented 
with numerical simulations of the electromagnetic wave 
transmission through the chain with a finite concentra- 
tion of impurities. The simulations demonstrated the 
development of a defect-induced polariton band within 
the forbidden polariton gap. In this paper, we present 
results of an analytical study of properties of this band. 
Neglecting spatial dispersion of phonons, we are able to 
obtain analytical expressions for DOS and the Lyapunov 
exponent for the band in the limit of large concentra- 



tion of impurities lo/ldef ^ I7 where Idef is the average 
distance between impurities, and Iq is the localization 
length of the single local polariton state. In this limit 
the method of microcanonical ensemble, complemented 
by expansion over (lo/ldef) ^ , gives the DOS of an ef- 
fective system with a uniform distribution of impurities, 
and describes its renornialization caused by local fluctu- 
ations of impurity positions. It also reveals localization 
of the states of the band, which are all localized due to 
the one-dimensional nature of the model. The properties 
of this polariton band are found to be drastically diffcp-, 
ent from that of the band produced by two-level atoms.tZl 
The analytical results show excellent agreement with nu- 
merical calculations and lead to suggestions regarding ex- 
perimental observation of the impurity induced polariton 
bands. 



II. THE MODEL OF ONE-DIMENSIONAL 
POLARITONS AND THE METHOD OF 
CALCULATIONS 

1. The model 

The main objective of the present paper is to study 
properties of the impurity induced polariton band aris- 
ing within the polariton gap of a one dimensional chain 
of atoms with dipole moments P„, where n refers to a 
position of an atom in the chain. Dynamics of the chain 
interacting with a scalar "electromagnetic field" E{x) is 
described within the tight-binding approximation with 
an interaction between nearest neighbors only, 

(nl - + + ^ aE{xn). (1) 

where $ is the parameter of the interaction between the 
dipoles, a is a coupling parameter between the dipoles 
and the electromagnetic field, and fi^ represents a site 
energy. We shall assume that the chain is composed of 
atoms of two different types, which differ in their site 
energy only. In this case fl^ can be presented in the 
following form: 

nl = nl + Afi^cn, (2) 

where Afi^ = — JIq, and flf are site energies of 
the respective atoms, and c„ is a random variable which 
takes values and 1 with probabilities I — p and p, re- 
spectively. We assume that p ^ 1 so that ilo is at- 
tributed to host atoms, whereas fli corresponds to the 
impurities. Eq. (|^) can be interpreted in terms of both 
exciton-like and phonon-like excitations. In the latter 
case, AfJ^ = {l-Mijnp/Mhost)i^'^, where Mhost and M^^p 
are the masses of host atoms and impurities, respectively. 

Polaritons in the system arise as collective excitations 
of dipoles (polarization waves) coupled to the scalar elec- 
tromagnetic wave, E(xn). The equation of motion of the 
electromagnetic component is given by the Maxwell equa- 
tion: 
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E{x) 



= — 47r— Pn6{na — x), 



(3) 



where c is the speed of the electromagnetic wave in vac- 
uum and a is the distance between two nearest-neighbor 
dipoles. Eqs. and present a microscopic descrip- 
tion of the transverse electromagnetic waves propagating 
along the chain in the sense that it does not make use of 
the concept of the dielectric permeability, and takes into 
account all modes of the field including those with wave 
numbers outside of the first Brillouin band. 

Considering propagation of the electromagnetic wave 
as a sequence of scattering events separated by free prop- 
agation, one can present Eq. (^) in a discrete form. The 
entire system of dynamic equations (|^) and (|3|) can then 
be transformed into a transfer-matrix form: 



Vn+l 



(4) 



Vectors u„ with components (P„, Pn+i, E„, E!^/k), where 
En and E'^ are electromagnetic field and its derivative at 
the n-th site, represent the state of the system. In gen- 
eral, the matrix ■?„, which describes the jEekLpropagation 
between adjacent sites, is a 4 x 4 matrixE3"E3 This makes 
an analytical consideration of Eq. very cumbersome, 
even in the case of a single impurity. At the same time, 
numerical simulations carried out in Refs. showed 
that most of the features of the impurity induced polari- 
ton band are well reproduced even when the interaction 
between the dipoles (spatial dispersion) is neglected. In 
this approximation, the two first components of the dy- 
namic vector Vn are easily expressed in terms of the other 
two, and the transfer matrix fvt becomes a 2 x 2 matrix 
of the following form: 



cosfca sinfca 
- sinfca -|- (5n(ka) cos ka cos ka + I3n(ka) sin ka 

(5) 

The parameter /?„ 



I3n 



(6) 



with (P = Ana/ a represents the polarizability of the n-th 
dipole. 

The main object of interest in studies of one- 
dimensional disordered systems is the Lyapuno'v— expo- 
nent (LE), A, which in our case can be defined ag23 



A = lim — In 



111 



\vo\\ 



(7) 



where L is the total length of the chain consisting of 
N atoms. It is well-known that the quantity defined 
by Eq. ^ is non-random (self-averages). It character- 
izes the spatial extent of the envelope of system eigen- 
states, all of which are localized in one dimension (see 



Ref. |2^ and references therein). The same quantity also 
describes the typical value of the transmission coefficient, 
T, of an external excitation incident upon the systemi 
Ttyp ~ exp(— AL), and through the Thouless relationcfl 
it determines the DOS of the system. Calculation of LE 
in the spectral region of the polariton gap of the pure 
system is the main task of this section. 



2. The Method of Microcanonical Ensemble 

The major problem with calculating of LE or DOS in 
the region of impurity bands is that a simple concentra- 
tion expansion is not able to describe the impurity band. 
Electron and phonon impurity bands have been inten- 
sively studied in the past (see, for example, Ref. In 
the limit of extremely small concentration of impurities, 
(plo) I a <C 1, it was proved possible to provide a regular 
systematic of the arising stateg|-,and use statistical ar- 
guments to describe their DOS.eI Another method used 
for this type of calculations employed the so-called phase 
formalism (sec Ref. ^ and references therein) . This ap- 
proach allowed one to calculate DOS of impurity bands 
in different spectral regions, including exponential be- 
havior at the tails of the band for both small and large, 
(pZo)/a ^ 1, concentrations. In the case of local po- 
laritons the localization length, ^q, of an individual local 
state can be so large that one has to deal with the situa- 
tion described by the latter inequality even in the case of 
a rather small concentration of impurities. Therefore, 
in this paper we focus mainly upon the properties of 
the well-developed polariton impurity band when indi- 
vidual local states are overlapped. To this end we use 
the method of the microcanonical ensemble, which was 
first suggested for analytical calculations of LE in a onfe 
dimensional single-band model of a disordered alloy.cS 
The advantage of this method over other methods of DOS 
calculations is that it allows one to calculate simultane- 
ously both DOS and the localization length. Its main 
shortcoming is that, as we shall see later, this method 
is a version of "effective media" methods, and as such 
it is unable to describe DOS near fluctuation boundaries 
of the spectrum. At the same time the results obtained 
allow a clear physical interpretation and, as follows from 
comparison with numerical simulations, they quite ac- 
curately describe properties of the bulk of the impurity 
band. 

The starting point of our calculations is the general 
definition of LE given by Eq. (0). In spite of LE be- 
ing a self-averaging quantity, it is convenient in practical 
calculations to perform averaging over the random con- 
figurations of the impurities. The regular ensemble of 
the realization is described by the fixed concentration 
of the impurities p, while their total number can vary 
from realization to realization. Therefore one can distin- 
guish two causes for fluctuations: (1) local arrangements 
of the impurities and (2) the total number of the impu- 
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rities. The main idea of the method of microcanonical 
ensemble is to reduce the finite size fluctuations in the 
system, ehminating the fluctuations of the total number 
of impurities. Such an ensemble with a flxed number 
of impurities is called the microcanonical ensemble in 
analogy with statistical physics. At the same time the 
result of the averaging in the limit L — > cx) should not 
depend upon the type of ensemble used by virtue of the 
self-averaging nature of LE. The key idea of the micro- 
canonical method is based on the assumption that with 
one cause of fluctuations eliminated, one can obtain re- 
liable results when the microcanonical ensemble average 
of (ln(...)) is replaced by ln(...). Such a substitution gives 
an exact result in the case of commutating matrices, and 
leads to an excellent agreement between analytical calcu- 
lations and simulations in. the case of 2 x 2 matrices with 
a single band spectrum.t^ 

Using the microcanonical ensemble, one can evaluate 
the average over all matrix sequences using the following 



expression derived in Rcf. 23 



/ N 

n 



{pN, 



gpN 



ipN)l dxP^ 



(fo -I- XTl) 



N 



(8) 



where tq and fi are host and defect matrices, respectively 
and a: is a free parameter. The derivative in Eq. then 
can be presented in the form of the Cauchy integral 



/ N 

n 



N 



pNj 2mJc xP^+^ ^' 



(9) 



where the contour of integration C is taken on the com- 
plex plane around x = 0, ly (x) is the largest eigenvalue 
of the matrix fo -I- xti, and the matrix D(x), 



D{x)^v-^{x) (fo+xfi) 



N 



(10) 



has eigenvalues not exceeding 1 in absolute value. In the 
limit of large N the above integral can be evaluated by 
the saddle point method. As the result we arrive at the 
following expression for the complex valued LE, A: 



A = lim — In 



largest eigenvalue ^ JJ^ f „ 



In //(xo) — (1 — p) In 



l-p 




where xq is defined by the saddle point equation 
dhiv {x) 



dhix 



p. 



(11) 



(12) 



The real part of A given by Eq. ( [Til ) represents LE, while 
its imaginary part according to ThoulesaHJ gives the in- 
tegral density of states N{uj) in the impure system: 



A = Re 



N{lo) 



1 



Im 



(13) 
(14) 



The eigenvalues of the matrix tq -f xti can be found 
from the equation 



l + x 



1 + x 



K (x) — 2 cos (fca) 



l + x 



(15) 



where (3^ and [3i are polarizabilities of the host and defect 
atoms, respectively. 

It is convenient to rewrite Eqs. ( pT| ) and ( p^ ) in terms 
of a new variable y — x/ [1 + x): 



«(y) = (1 - vy{v) + 



1 



(1 - y)v{y) 



Ko{(-o) — 2cos(fca) -|- (i^fcasin(fca)- 



F{lo) — d^kasm{ka) 



(Cj2 - %) (a,2 - J7f ) 



ko(cc;)+j/F(l^), (16) 
(17) 
(18) 



"o 



In the long-wave limit e = fca ^ 1 functions F{tj}) and 
Ko(ijj) can be simplified and presented in the form, which 
clarifies their physical meaning: 

Ko{uj) ~ 2 4- 6^7(0;), 
where /(w) is defined according to 



/(^) 



(19) 



and represents the difference /3i — f3o between polarizabil- 
ities (|) of the impurities and host atoms. Function j{oj) 
defined as 



7(tj) = 1 + 



(20) 



is the longwave dielectric function of the pure chain. 
Similarly, Eq. (|l^) transforms into 

Xip) = \ii[il-yo{p))iy{yoip))] 

pln^ + (l-p)lnl-^"(P) 



p l-p 
with the saddle point equation reading as 



± 



y(i - y)F 



^J(no + yFY 



(21) 



(22) 



The choice of the sign in Eq. ( p^ ) is determined by the 
requirement to use the greatest of the eigenvalues ly. In- 
troduction of the new variable y turns k (y) in a linear 
function of y essentially simplifying future calculations. 
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III. THE LYAPUNOV EXPONENT AND THE 
DENSITY OF STATES 

3. Boundaries of the impurity band 

Dielectric function "f{uj) determines the frequency re- 
gion of the polariton gap: for < oj^ < fl^ + d? it is 
negative, and hence, propagating modes do not exist in 
this region. We shall assume that the defect frequency ^\ 
obeys the inequality fig < 17^ < JIq+c?^, so that the impu- 
rity induced band develops inside the gap of the original 
spectrum. As we already mentioned, the approach we use 
in the paper belongs to the class of effective-medium ap- 
proximations since we neglect here certain kinds of fluc- 
tuations in the system. It is natural to expect, therefore, 
that within this approach the impurity band would have 
well-defined spectral boundaries outside of which the dif- 
ferential DOS remains exact zero. Our first goal is to 
determine these boundaries and find the concentration 
dependence of the width of the impurity band. The dif- 
ferential DOS, p{lu) = dN{u:)/dcu, takes on non-zero val- 
ues when X{uj) acquires a non-constant imaginary part. 
Rewriting Eq. (E^) in the form 



Xiu) = i In 
a 



yo 



yoF± 



yo(i - ya)F 



pin — + (1 ~ p) In 
p 1 



P - yo 

1-2/0 



P 



(23) 



one can see that in the case of real yo, ImA(w) can be 
either zero or tt (the latter happens when the argument 
of In is negative). In both cases differential DOS is ob- 
viously zero, and it can only take on a non-zero value 
when yo becomes complex. Eq. (^^, which defines yo, is 
formed by a polynomial of the third order with real coef- 
ficients; therefore it has either three real roots or one real 
solution and a complex conjugated pair. In order to de- 
scribe formation of the polariton band, one has to select 
the root, which has a complex component in a certain 
frequency interval and yields positive LE. Comparison 
with numerical simulations shows that the choice of only 
one of three saddle point, according to above mentioned 
criteria, produces the correct description of the impurity 
band. It is not difficult to show that at the frequency 
where yoi^) becomes complex, the derivative dyo/duj di- 
verges. This fact gives us an explicit equation for the 
spectrum boundaries. From Eq. ( |2^ ) one can find that 
the divergence occurs when yo satisfies the equation 



ylF [F{1 -p) + Ko] + yl [koF{1 - 3p) + 



4] 



+yoHo {F - 2nl + 4) + p{kI - 4) = 0. (24) 

Eqs. (^2|) and ( p4[ ) define the concentration dependence 
of the boundaries of the polariton impurity band. An 
approximate solution of this equation can be obtained as 
a formal series in powers of parameter e = ka. It will be 
seen, however, from the results that the actual expansion 



parameter in this case is \/plo *C 1. The solution of Eq. 
(E2) with the accuracy to can be obtained in the form: 



y=p-p{\-p) 



+ p{l-p) 



f 



2(p/-7)/(l-2p)-/p(l-p) 
8(p/-7)' 



(25) 



where y — p is the only non- vanishing zero order approx- 
imation for y. Since two other solutions,?/ — 0, corre- 
spond to a singular point of the integral (0) , the chosen 
solution y = p represents the only saddle point accessible 
within our perturbation scheme. Using additional crite- 
ria outlined above we verify that the solution given by 
Eq. ( p5[ ) correctly reproduces the behavior of LE. Sub- 
stituting Eq. ( p^ ) into Eq. (|2^) one finds the complex 
LE in the long wavelength approximation: 



C 



Vpfi^ 



> 8[p/(c.) - 7(^)] 



(26) 



According to Eq. (|26[), \{uj) acquires an imaginary part 



at frequencies obeying the inequality 
p/(c^)-7(c^) <0, 



(27) 



The boundaries are determined by the respective equa- 
tion pf(uj) — ^{uj) — 0, which coincides with the long 
wavelength limit of Eq. (p^. This equation determines 
two points where pf{uj) — 'y{i-^) changes sign to negative: 



2 



-inl + nl + d^) 



(28) 



d^-niY 



The first of these solutions belongs to the initial polari- 
ton band-gap and as such represents the low-frequency 
boundary of the new impurity band. The second one lies 
outside of the gap and is a bottom frequency, modified 
by impurities, of the upper polariton branch of the ini- 
tial spectrum. These two frequencies, however, are not 
the only points where the expression pf{uj) — ^{lo) turns 
negative. Two other points are 



OJ., 



(29) 



and the change of the sign at these points occurs through 
infinity of pf{i-o) — ^{lo) rather than through zero. These 
two frequencies do not depend upon the concentration of 
impurities and present, therefore, stable genuine bound- 
aries of the spectrum. This property of rJo and lli is due 
to their resonance nature (they correspond to the poles 
of the respective polarizabilities) and disappears when, 
for example, the spatial dispersion is taken into account. 
At the same time, the numerical simulations of Ref. |2|, 
in which the spatial dispersion was taken into considera- 
tion, indicate that the shift of these frequency from their 
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initial values is negligibly small for realistic values of the 
inter- atom interaction parameter even for a relatively 
large concentration of impurities. 

These four frequencies set the modified boundaries of 
the initial polariton spectrum of the pure system (cj^;, 
and boundaries of the newly formed impurity band 
{wfy^, cofi). The lower boundary of the forbidden gap of 
the crystal, fig, is not affected by the impurities - the 
singularity in the polariton DOS at this point survives 
for any concentrations of defects. The upper bound- 
ary of the band gap, LUpu, shifts toward higher frequen- 
cies with the concentration and when p = 1, it coin- 
cides wit h the up per band boundary of the new crystal, 
= + (P. Frequencies ujn and uJiu give approxi- 

mate values for the lower and upper boundaries of the 
impurity induced pass band which arises inside the orig- 
inal forbidden gap fio < < ^iu < + d'^- The 
impurity band grows asymmetrically with concentration: 
while its lower boundary moves towards f2o with an in- 
crease of the concentration, the upper edge remains fixed 
at LUiu — ^li- Such a behavior of the impurity band agrees 
well with our numerical results. c3 

The width of the impurity induced band defined in 
terms of squared frequencies 
from Eqs. and @ as 



-UJZ 



-ujfi can be found 



1 



(U 



(30) 



For Oi not very close to the upper boundary of the initial 
polariton gap, J^l = y^l^ + d^, the linear in p approxi- 
mation is sufficient to describe the concentration depen- 
dence of the band width 



A2 



d^ 



(31) 



When rii, however, happens to be close to JIl, a crossover 
is possible from the linear dependence (|3^) to the square 
root behavior: 



A^^^^dHnl-nl)p. 

The condition for such crossover to occur is 

(nl + d^-nl)^ 



(32) 



Ad^nj-ni) 



4- The Lyapunov exponent and the density of states far from 
the spectrum boundaries 



Xito) 



P(l -P)/^ ^ 
" 8 (p/ - 7) c 



e + 0(e'), 



(33) 



1 U} 



TT C 

for forbidden bands, (f2oi'^ii)i i^ii'-^pu), we obtain 



X{uj) = ^ jVp./ - 7 - 
N{uj) = 0. 



8 {v! - 7) 



(34) 



One can see from Eqs. (^3|) and (|3J) that within allowed 
bands, DOS appears in the zero order of the formal ex- 
pansion parameter e, while LE starts from the first order. 
For the forbidden bands the situation is reversed: LE 
contains a term of zero order in e, while DOS in this or- 
der disappears. This observation suggests a simple phys- 
ical interpretation of the microcanonical approximation 
in conjunction with the expansion over e. Analysis of 
the series in this parameter shows that the actual small 
parameter is a/{plo) = Idef/^o, where Idef is a mean dis- 
tance between impurities. The zero order expansion in 
this parameter can be interpreted as a uniform contin- 
uous {Idef 0) distribution of impurities with concen- 
tration p. The results in this approximation then cor- 
responds to the polariton impurity band which would 
exist in such a uniform system. The parameter Idef /^o in 
this case is a measure of disorder in the distribution of 
impurities, which leads to the localization of excitations 
described by the DOS in Eq. (|3^) with the localization 
length, l^^ — A, presented in the first line of the same 
equation. 

Let us now consider frequencies in the interval uj G 
(0, f^o) U (il^,oo), where pass bands of the pure chains 
composed of either host or impurity atoms overlap. The 
DOS in this region can be written in a physically trans- 
parent form 



N(w) - ^{\-p)Nl{^)+pNl{w)- 



(35) 



where Nq{uj) and Ni[uj) are integral DOS in pure chains 
containing host atoms or impurities, respectively. In the 
remaining portion of the initial spectrum lu £ (iOup, ^l) 
(but not very close to the boundary tOup, where the ex- 
pansion ceases to be valid) DOS can be presented as 
N{lu) = ^(1 - p) Nq{uj) — pl\{uj), where is the pen- 
etration length through the polariton gap of the 100% 
impure chain. 

Our main goal, however, is obtaining DOS and LE 
of the impurity induced band, oj S {oJii,^i). Near the 
center of this region, the DOS can be presented in the 
following form 



Using Eqs. (|T3|), (|T|) and ^ one can calculate LE 
and the integral DOS for different regions of the spec- 
trum. For allowed bands, (0, fio); {'^Ut^i), (wp„,oo), 
one has 



TT/q 



4|7|^c 
pd"^ 



O 



(36) 
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where ujc is the center of the impurity band, which in the 
hnear in p approximation is 



(37) 



where A^,^ is the width of the band in terms of squared 
frequencies defined in Eq. (|l|). The first ter m in Eq. 
( Pq ) represents the total number of states between uJa 
and the center of the band; it is interesting to note that 
this number does not depend upon concentration of im- 
purities. The coefficient at the second term gives the 
differential DOS at the center and can be rewritten as 



4|7h 



ttIqpcP ttIqS ' 



(38) 



where 6 = Af^/2aJc is an approximate expression for 
the impurity band width (5 ~ f2i — lou. This DOS has 
a simple meaning of the average density of states uni- 
formly distributed through the entire band over the dis- 
tance equal to lo/2. In one dimensional uniform systems 
the wave number of the respective excitations, k, is sim- 
ply connected to N (u): k — ttN{uj). The differential 
DOS would in this case be proportional to the inverse 
group velocity. Accordingly, l/np^LJc) given by Eq. (^) 
can also be viewed as a group velocity, of excitations 
in the center of our impurity band in the case of the uni- 
form distribution of impurities. The expression for v can 
also be presented as 



3/2 



pc < C, 



(39) 



which demonstrates that polariton excitations of the im- 
purity band have much slower velocities not only com- 
pared to c but also to the velocities at both regular po- 
lariton branches. 

Expanding Eq. ( p3| ) for LE about the center of the 
band, we obtain a parabolic frequency dependence of the 
localization length of the impurity polaritons 



1{lo) = X{uj)-^ ^2loiplo) 



1 - 



2O72 



? ? \ 2 



<p 



(40) 



One can see from this expression that / (w) reaches its 
maximum value 2pl'^ at the center of the band. It is im- 
portant to note that the localization length here grows 
linearly with the concentration, whereas it is LE that 
grows with concentration for frequencies outside of the 
impurity band. An increase of the concentration of the 
impurities also results in fast (cx 1/p^) flattening of the 
maximum in the localization lengtb—the fact we first no- 
ticed in our numerical simulations E3 

The frequency dependence of LE, defined by Eqs. ( |2l| ) 
and (p^, is shown in Fig. 1. For frequencies correspond- 
ing to the impurity band, LE drops sharply (the local- 
ization length increases); it then diverges at the upper 



boundary of the impurity band. This divergence has the 
same origin as LE divergence at the lower boundary of the 
band gap of the original crystal, where due to the speci- 
ficity of the polariton spectrum, the wave vector becomes 
infinite. The concentration dependence of LE and inte- 
gral DOS for some frequency cjo G (f^o, f^i) is shown in 
Fig. 2. For p = this frequency belongs to the forbidden 
gap, thus DOS is zero. With an increase of the concen- 
tration LE decreases (the localization length increases). 
For the concentration when the lower boundary of the 
impurity band crosses ujq, DOS becomes non-zero and in 
\{uj) the crossover between behaviors described by Eqs. 
( p^ ) and ( p^ ) occurs. In Fig. 3 we compare the results of 
our analytical calculations with numerical simulations of 
Ref. ^ The comparison shows an excellent agreement 
between numerical and analytical results, confirming the 
validity of the microcanonical method in the considered 
limit. In addition, since numerical results were obtained 
with the spatial dispersion taken into account, the com- 
parison shows that the model with dispersionless phonons 
produces reliable results for LE even not very far away 
from the spectrum boundaries. 



5. The solution in the vicinity of the spectrum boundary. 
N onanalytical behavior. 

The results obtained in the previous section are clearly 
not valid for frequencies close to the band boundaries oju 
and ojpm where pf — j = and the second term in the 
expansion ( |25| ) diverges. It is well known that perturba- 
tion expansions in disordered systems usually fail in the 
vicinity of boundaries of the initial spectrum of the sys- 
tem unperturbed by disorder .Ej The regions in the vicin- 
ity of these special frequencies uu and Upu require special 
consideration. Attempting to regularize our e-expansion, 
we shall seek corrections to the zero order solution in the 
form 



y^p + Be" 



(41) 



admitting a possibility of fractional values of a, and, 
hence, non-analytical behavior of the solution. We also 
introduce a new variable ( which determines the proxim- 
ity to either of two frequencies uja or iUpu' 



(42) 



Substituting these expressions into Eq. ( [2^ ) and equat- 
ing the lowest order terms we see that the equation can 
only be satisfied for a = a' = 2/3. In this case we find 
that the parameter B introduced in Eq. (^l]) obeys the 
equation 



2VCTBT 



B = 0. 



(43) 



The substitution of y given by Eq. (^Tj) in the condition 
for the spectral boundaries presented by Eq. (|3) allows 
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one to obtain an exact expression for the renormalized 
boundary 



Pf -1 



2/3 



4/3- 



^2/3 



(44) 



that modifies Eq. (|2^). The shift is small by virtue of 
smallness of e^/^ . In the lowest order in e, the new posi- 
tions of the band edges for small p are 



= n\ 



^pu L 



2/3 



(45) 



1 - 3 



16plo 



1/3 



where we again encounter a/{plo) as a true small parame- 
ter of the expansion. It is interesting to note the different 
character of non-analytical corrections to the positions of 
the boundary of the impurity band ujfi and the bottom 
of the upper polariton branch Wp^. They have fractional 
concentration dependence with the correction to be- 
ing much stronger in the limit a/{plo) «C 1. 

Now let us consider the modification of LE and the in- 
tegral DOS due to the nonanalyticity. Substituting Eq. 
( ^ ) into Eq. (|2^) the complex LE can be written as 



A: 



P(l -P)/ 4/3 

2B ^ 



{pf - 7) + BfeV^ 



(46) 



This expression explicitly demonstrates the crossover be- 
tween analytical and nonanalytiacal behavior. When fre- 
quency uj lies far from the initial boundary \pf — 7I ^ 
Bfe^^^, one recovers the term proportional to e with the 
same coefficient as in Eq. (p6[), whereas in the opposite 
limit \pf — 7I ^ A/e^/-^, when we approach the bound- 
ary, the leading term gains fractional power and becomes 



oc e 



4/3 



At the vicinity of the renormalized spectral boundary, 
Eq. (^, DOS p{uj) = to the left of uj^i and for u; > uju 
it can be obtained as 



p{uj) 



3 



1/2 



1 



-2^,1/2' 



(47) 



where 7 and Iq are evaluated at w = Oi, and we have ne- 
glected the renormalization of the boundary when calcu- 
lating the coefficient in p (w) . In this approximation the 
frequency and concentration dependence of DOS does not 
change as compared to the one obtained from the non- 
renormalized expansion (|3^). It is interesting, however, 
that the renormalization brings about an additional nu- 
merical factor of 3 in Eq. (^^, which is absent in the non- 
renormalized expression. This frequency dependence is 
typical for excitations with the quadratic dispersion law 
in the long-wave approximation. The main characteris- 
tics of this dispersion law is the effective mass to, which 
can be found from Eq. (47) as 



97 _ 972 
2pdHl ^ 2^ 



(48) 



It is interesting to compare this expression with the ef- 
fective mass of the upper polariton branch of the pure 
system at the bottom of the band 2too = {ojL/dcf'. The 
two expressions have a similar structure if one introduces 
a "renormalized speed of light" 



c ~ c\/pl^y^. 



(49) 



Though the introduced parameter c does not have a di- 
rect meaning of the speed of the excitations, it shows 
again that the excitations in the impurity polariton band 
are considerably slower than their regular counterparts, 
with a similar dispersion law at the spectrum boundary 
as well as at the center of the band. Unlike, however, the 
center-of-band situation ( p9| ) the renormalized velocity 
at the edge is proportional to the square root of concen- 
tration. 

LE in the vicinity of Zjh is represented by different ex- 
pressions for frequencies below and above ujn respectively 



2lo \Aplo 



1 



2/0 V4p'' 



1 



1/3 



1/3 



12 (7^0)'^' 
d (4Zop)'/' 

12 7^0 



(50) 



rf2 



(4/op) 



1/3 



refiecting a discontinuity of its frequency derivative at the 

spectrum boundary. LE itself is, of course, continuous at 

1/3 

LUii, giving rise to the localization length I = 21^ (4pZo) 

at the edge of the band. It is interesting to compare 
this expression with the localization length at the cen- 
ter of the impurity band, Eq. (|40|). They both grow 
with the concentration but the latter one is much smaller 
and demonstrates slower fractional concentration depen- 
dence. 

At the upper impurity band edge fli integral DOS, 
N{uj), diverges causing much stronger singularity in the 
differential DOS than at the lower boundary luh 



p{uj) = 



dujp^/^ 



TTC — LlP- 



3/2 ' 



(51) 



which is typical for DOS in the vicinity of resonance fre- 
quencies. Comparing this expression with a similar for- 
mula for a pure chain, one can again interpret this result 
as a renormalization of the velocity parameter c by the 
concentration of impurities, which is different from that 
presented by Eq. (^9|) by a numerical coefficient only. 

The last spectrum boundary is the bottom of the up- 
per polariton branch O^. The spectrum in the vicinity 
of this frequency exists in the absence of the impurities, 
which are responsible for two effects in this region. First 
they move the boundary from f2i to higher frequencies 
[Eq. (IH)], and second, they increase the effective mass 
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of the upper polariton branch, such that the differential 
DOS in the frequency region above cD2„ becomes 



p{lu) 



1 



up 



(p)) 



1/2 



(52) 



This expression can also be interpreted as a renormaliza- 
tion of the speed of light c. 



IV. CONCLUSION 

In the paper we have presented a detailed study of 
an impurity induced band of excitations, which arise in 
the gap between lower and upper polariton branches of 
a linear chain with dipole active atoms. We have also 
studied impurity-induced effects on properties of the reg- 
ular polariton branches. The method of microcanonical 
ensemble in conjunction with the expansion in the pa- 
rameter Idef/io 1) where Idef is the average distance 
between the impurities and is the localization radius 
of a single local polariton state, produces a clear phys- 
ical description of the excitations, and shows excellent 
agreement with the results of numerical simulations. In 
the zero order of this expansion, we recover DOS and 
the dispersion law of excitations in the system with uni- 
form continuously distributed impurities. Corrections to 
this solution describe effects due to local fluctuations in 
the positions of impurities such as a finite localization 
length of the excitations, a renormalization of the spec- 
tral boundaries and the effective mass of the excitations. 
The parameter Idef/lo, therefore can be considered as a 
measure of disorder in the system. 

The main result of the calculations is the demonstra- 
tion of the existence of the polariton band formed by 
impurities without inner optical activity. Justification 
for naming these excitations as polaritons lies in the fact 
that they carry an electromagnetic field coupled with me- 
chanical excitations of the chain. The dispersion law of 
the excitations at the lower frequency spectral boundary 
resembles that of the upper regular polariton band, but 
with a significantly (by a factor of p^^^ <€. 1) reduced ef- 
fective mass. At the higher frequency edge of the band, 
the wave number diverges in a manner similar to the 
regular lower polariton branch with a velocity parameter 
again reduced by the same factor. The group velocity 
of the excitations near the center of the band, as well as 
is the localization length, is proportional to p. The lat- 
ter, however, demonstrates a parabolic frequency depen- 
dence with the curvature of the parabola falling off with 
increase of the concentration as Excitations consid- 

ered in this paper are drastically different from impurity- 
induced polaritons studied in Ref . The authors of the 
latter paper considered excitations of an ordered chain of 



two-level atoms embedded in a polar 1-d crystal. Most 
significantly, the effective mass of polaritons of Ref. |l^ 
is negative at the long-wave boundary of the spectrum 
while our excitations have a positive effective mass in 
this region; concentration dependencies of the effective 
mass and the band width also differ significantly. 

The regular expansion in powers of the parameter 
Idef/lo produces diverging expressions at the boundaries 
of the zero order spectrum. Allowing for fractional pow- 
ers in the expansion, we obtained a finite regularized 
expression for the density of states and the localization 
length at the boundaries. It is interesting to note that 
renormalized DOS at the lower boundary of the impu- 
rity band differs from its initial form only by a numerical 
factor of 3, while the position of the boundary is shifted 
toward lower frequencies. 

Experimental significance of the considered excitations 
is affected by two factors: absorption due to different 
kind of anharmonic processes, and the one-dimensional 
nature of the considered model. It is generally accepted 
(see for instance Ref. ||) that one-dimensional models give 
a fair description of tunneling in the limit of small con- 
centrations {Idef/lo ^ !)■ We deal with the opposite 
limit; however, our zero order results could describe the 
dispersion law of excitations in the real three dimensional 
medium propagating in one specified direction (for exam- 
ple, the direction with the highest symmetry). Disorder 
in this case would lead, of course, to scattering and devi- 
ation from one-dimensional geometry, but effects due to 
disorder are much weaker in three dimensions, and would 
not probably inhibit an observation of a defect induced 
transparency in the frequency region of a polariton gap. 

One can give a simple estimate for the role of absorp- 
tion comparing the absorption coefficient with the width 
of the impurity induced band. If the latter is greater the 
band would manifest itself in transmission (or reflection) 
experiments. Using Eq. ( |3l| ) one can obtain the following 
estimate for the admissible absorption coefficient a: 



a < 



(P{nl-nl)p 



if the defect frequency fii is not very close to the upper 
boundary of the gap fl^ or 

in the opposite case. The latter situation is much more 
favorable because of the square root concentration de- 
pendence. For example, for frequencies of the order of 
1000 cra~^ and volume concentration of 0.1%, which 
correspond to linear concentration entering our formu- 
las p ~ 10%, one has from the last expression a <C 300 
CTO~^, while the typical value of absorption for, e.g., al- 
kali halides is ~ 100 cm~^. This estimate shows that an 
observation of the transmittance through polariton gaps 
due to impurity polaritons is feasible. 
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FIG. 1. LE for a chain with a defect concentration of 10 % 
(solid line) in comparison to a pure system (dashed line). 



FIG. 2. Dependences of LE (sohd line) and DOS (dashed 
line) on concentration for a frequency in the interval 
f2o < w < fii. 



FIG. 3. Comparison of LE calculated with (dashed line) 
and without (solid line) spatial dispersion in the vicinity of 
Qi. The concentration of defects for both curves is 1 %. 
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